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1.  INTRODUCTION 

Several  important  mechanisms  of  air-ocean  inter¬ 
action  for  the  El  Nino  and  Southern  Oscillation  (FNSOl 
phenomenon  have  been  developed  in  the  past  dev  tskr 
ocean  wave  propagation,  deby -oscillator.  two  equilib¬ 
rium  states,  and  air-ocean  coupled  instabilities  Due  to 
crude  parameterization  of  thermodynamical  processes  in 
bosh  atmosphere  and  ocean.  these  theories  either  cannot 
explain  the  transition  between  El  Nino  and  La  Nina 
(e  g.,  coupled  instability  theories),  or  need  an  artificially 
setting-up  criterion  to  make  such  a  transition  (c.g..  slowly 
propagating  oceanic  Rossby  wave  theory).  The  ineg- 
ulnrity  of  ENSO  implies  that  the  ENSO  events  cannot 
he  explained  as  n  pure  wave  propagation.  More  detailed 
research  on  thermodynamics  in  both  ocean  and  atmos¬ 
phere  is  needed  before  we  run  the  oceanic  and  atmo¬ 
spheric  GCMs. 

2.  RELATIONSHIP  BETWEEN  PERTURBATIONS 
OF  SSI  ANDOMI  DEPTH 


The  thermodynamics  ol  upper  ocean,  commonly 
used  in  the  coupled  models,  can  be  summarized  bv 
(Hirst,  19X6) 

(?r.  r>?.  dT. 
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where  r u'„  A'.  arc  perturbations  of  SST.  zonal  cur¬ 
rents.  and  OMI.  thickness.  PTJcx  is  the  mean  zonal  SS I 
gradient,  k  and  a  arc  two  parameters.  If  k  is  large. 
Eqi  1 1  becomes 

T,  -  eA',  (21 

The  expressions  show  that  the  mcrease  of  the  OML 
thickness  (A-. >  0)  leads  to  the  increase  of  SST  ( 
dTJft  >  0.  or  r.  >  0)  Is  this  type  of 

thermodynamics  correct  ?  Not  really.  It  is  correct  onlv 
for  a  cooling  OML.  where  the  outgoing  heat  flux  is 
greater  than  imraning  heat  flux.  The  thicker  the  layer, 
the  less  cooler  of  the  layer.  Thus,  the  relationship  be¬ 


tween  perturbations  of  SST  and  OMI.  thickness  ts  given 
t*y 
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For  a  warming  OML.  on  the  contrary,  the  thinner  the 
layer,  the  hotter  the  layer.  Therefore,  the  relationship 
between  perturbations  of  SST  and  OMI  thickness  is 
written  by 


an 

>v 


-  /.'  - 


(.Vi) 


The  fact  >hat  only  (3a>  is  chosen  is  r  therm-dynamical 
component  for  the  ocean  pan  in  coupled  air -ocean  mod¬ 
els  makes  the  current  ENSO  theories  quite  incomplete 


X.  AN  OML  SWITCHER 


Arguments  are  cast  in  terms  of  simple  mixed  layer 
models  (Chu  cl  aL.  1990;  and  Chu  and  Garwood.  I^l), 
where  it  s  assumed  that  the  temperature  and  velocity  are 
uniform  over  some  depth.  A.,  called  OMI  depth,  and 
that  the  penettation  depth  of  solar  radiation  is  much 
smaller  than  A..  With  these  assumptions,  one  can  write 
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where  T,  is  the  SST;  ft.  is  the  horizontal  velocity;  is 
the  characteristic  sea  water  density;  r  the  sea  water 
specific  heat  under  constant  pressure;  <J„  is  the  net  sur¬ 
face  heat  Bun,  downward  positive;  Q  ,  is  the  entrainment 
heat  flux  at  the  base  of  the  OMI  computed  by 


tM  -  - r  *>  ,s> 

where  k  and  T arc  the  entrainment  velocity  and  the 
lempciuture  at  the  base  of  the  OMI 

Entrainment  velocity  can  also  be  parameterized  in 
terms  of  OML  TKE  balance  If  salinity  effect  is  no 
considered  at  present,  sources  and  sinks  of  TKfc  at  ocean 
surface  arc  wind  work  (proportional  to  wind  speed 
cubed i  and  buoy  ant  damping  (oi  forsirR)  due  to  surface 
warming  (or  cooEng). 

tC,«i*  - 

H*  = - - -  (** » 
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where  C,  and  C*  are  tuning  cneffirirntx,  and  ts  the 
water  surface  frictiun  velocity,  which  h  proptutional  to 
the  air  surface  friction  velocity  u, 

A.  is  the  Heaviside  function  of 
K,u.'.  A.  K  in  OMI  sw  if  cbm 

whose  value  equals  I  for  the  entrainment  regime,  usually 
associated  with  strong  surface  wind  forcing;  and  cqunk 
0  for  the  shallowing  regime  (the  OML  thickness  taken 
the  Monin-Obukhnv  length  scale),  usually  associated 
with  weak  surface  wind  forcing 


4.  AlR-fX'EAN  FEEDBACK  MECHANISM 
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The  tropical  ocean  surface  genet  ally  receives  heat, 
i.e  .  0>  -  0  .  The  tropical  subcloud  layer  is  usually  domi¬ 
nated  by  the  mean  easterlies,  and  the  atmospheric  deep 
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convection  is  usually  located  .it  the  upward  branch  of  the 
zonal  circulation.  In  the  cast  <01  west)  of  the  deep 
consectivc  region.  the  surface  wind  speed  k  enhanced  <ot 
reduced),  ami  therefore  the  OMI  depth  is  increased  tor 
decrease*))  Under  weak  Bird  forcing,  the  decrease  of 
OML  depth  (V,  <.  Ol  leads  to  the  atu merit. ition  of  ^ST 
perturbation  [Crjt'i  -  0)  and  in  turn  t**  a  production  of 
the  3tm<*sphcfic  convection  to  the  wcM  of  the  deep- 
convectise  region  (Fig.  la).  However,  under  strong  wind 
forcing,  the  increase  of  OML  depth  (/)'.  >  0)  leads  to  the 
augmentation  of  *>ST  perturbation  (.  7*,/f!r > 0)  and  in 
turn  to  a  production  of  the  atmospheric  convection  to  the 
cast  of  the  dccp-convcclivc  region  (I'ig. lb). 


Fie.  I.  Vimospheric  onscctmn  generated  by  air-ocean 
feedback  under  t a)  weak  wind  forcing,  and  (b)  the  strong 
wind  forcing. 


the  AOS  HE  model  pmltns  the  generation  of  unstable 
eastward  propagating  mode  with  phase  speed  (03-1.5 
m  s)  and  the  maximum  growth  rate  (I  S  yr)  appearing 
at  M  zonal  wavenumbers  (Fig.*).  where  the  rcawmb^ 
values  of  the  model  parameters  for  the  entrainment  re¬ 
gime  of  OML  are  j  -  0  02. 
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Fig.  2.  Schematic  presentation  of  the  AOSHF  model. 
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5.  AOSHF  MODEL 

The  AOSIIL  mo<lcl  (Fig.2)  is  n  cnuplod  system  con¬ 
sisting  of  the  Wind-Induced  Surface  Heat  Exchange 
(WISHt)  model  (Yano  nnd  Emnnuel,  I99||  nnd  the 
Ocean  Mixed  Layer  (OML)  model  (Chu  and  Garwood, 
1991).  Above  the  air-ocean  interface,  it  reduces  to  the 
ordinary  WISHF.  model  (Yano  ind  Fmanucl.  |90|)  if 
the  SST  perturbation  vanishes.  Below  tne  air-ocenn 
interface,  it  is  the  ordinary  OML  model  (Chu  and 
Garwood.  I90l|  The  AOSHE  model  t»  solved  analyt¬ 
ically  as  an  eigenvalue  problem  (Chu.  |99|ab)  for  two 
different  surface  conditions:  strong  and  weak  surface 
wand  forcing.  For  strong  wind  forcine,  there  arc  two  im¬ 
portant  parameters:  y  -  .tir-occnn  surface  coupling  coef¬ 
ficient.  and  y,  —  entrainment  parameter  For  weak  wind 
forcing,  there  is  only  one  parameter,  y  When  the  surface 
wind  stress  is  strong. 
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i  l.  iinxiv  i  unucr  strong 
ind  lorcmg,  on  ,|1C  p)1rnmerer  j,r„w  ,h  rntc.  ,, ,  , 

nnu  (bj  pha.se  speed  </«,'  ')  for  ,  -  f).o;>. 


When  the  ocean  surface  is  under  weak  surface  wind 
forcing  (equivalent  to  strong  surface  warming), 

cM'tyf,dalPv.<s)*"iut  +  -  C2<*xhwOolPwacPv> <  0. 

the  AOSHE  model  predicts  the  generation  of  unstable 
westward  propagating  mode  with  phase  speed  (0.4  m/s) 
and  the  maximum  growth  rate  (3.5/yr)  appearing  at  the 
lowest  zonal  wavenumber  (Fig.4),  where  the  reasonable 
value  of  the  model  parameters  for  the  shallowing  regime 
of  OML  is:  y  -  0.03. 

This  implies  that  there  is  an  inherent  switcher  in  the 
OML.  Shifting  the  OML  from  one  regime  to  the  other 
regime  only  depends  on  the  model  variables:  ut}  v4,  /?„, 
etc. 
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Fig.  4.  Dependence  of  the  AOSHE  mode  (under  weak 
wind  forcing)  on  the  parameter  y:  (a)  growth  rate 
and  (b)  phase  speed  (ms  ')• 


ft.  EL  NINO/LA  NINA  CYCLE 

Based  on  the  AOSHE  model  results,  I  offer  a  new 
theory  for  the  La  Nina/El  Nino  transition  depicted  ns 
follows. 

(a)  La  Nina  Conditions  Prevailed  Stage 

Starting  from  the  typical  La  Nina  conditions,  atmo¬ 
spheric  deep  convection  (denoted  by  "Convection  A")  is 
developed  in  the  western  Pacific  (Fig. 3a).  The  surface 
winds  in  the  central  and  eastern  Pacific  are  enhanced 
due  to  the  same  direction  of  the  mean  easterlies  and  the 
easterlies  associated  with  the  zonal  circulation  induced 
by  Convection  A.  The  OML  is  under  strong  surface 
wind  forcing.  The  unstable  convective  disturbances  (de¬ 
noted  by  "Convection  B  )  are  generated  in  the  cast  of 
Convection  A  and  propagating  eastward  (Fig.3b)  with 
phase  speeds  0.5- 1.5  m/s,  and  maximum  growth  rate  (~ 
1.8  /year)  appearing  at  the  3-4  zonal  wavenumber. 

(b)  Mature  of  La  Nina 

As  Convection  B  continuously  propagates  (toward 
east),  grows,  and  reaches  the  eastern  Pacific,  the  La  Nina 
enters  its  mature  and  transition  stage  (Fig. 5c).  During 
this  stage,  the  strong  surface  wind  stress  is  over  the  con¬ 
tinent.  Convection  A  and  Convection  B  competes  each 
other  by  reducing  each  other's  low-level  convergence  and 
high-level  divergence.  The  possibility  of  disturbance 
generation  between  Convection  A  and  Convection  B  is 
greatly  reduced  due  to  this  competition.  Because  Con¬ 
vection  B  is  flowing  over  a  relatively  cool  ocean  surface 
(La  Nina)  with  not  too  large  growth  rate  ( -  1.8 /year), 

two  outcomes  arc  expected:  (I)  Convection  A  survives, 
Convection  B  weakens  and  disappear;  the  system  goes 
back  to  La  Nina  (Fig.5a).  (2)  Convection  B  survives. 
Convection  A  weakens  and  disappear;  El  Nino  will  take 
place  (Fig.5d).  This  uncertainty  makes  La  Nina  cycle 
very  irregular. 

(c)  Onset  of  El  Nino 

As  Convection  B  survives  (Fig. 3d),  the  surface 
westerlies  are  prevailing  over  the  tropical  Pacific.  They 
counterbalance  the  surface  mean  winds  (easterlies)  and 
make  the  total  surface  winds  very  weak.  This  will  shift 
the  OML  to  another  regime  (i.e.,  shallowing  regime). 
The  unstable  convective  disturbances  (denoted  by  "Con¬ 
vection  C")  arc  generated  in  the  west  of  Convection  B 
and  propagating  westward  (Fig.Se)  with  phase  speeds 
-  0.4  m/s,  and  maximum  growth  rate  (~  3.3/year)  ap¬ 
pearing  at  the  lowest  wavenumber.  Convection  C  is  a 
fast  growing  mode. 

(d)  Mature  of  El  Nino 

As  Convection  C  continuously  moves  toward  west, 
grows,  and  reaches  the  eastern  Pacific,  the  El  Nino  enters 
its  mature  stage  (Fig.5f).  During  this  stage.  Convection 
C  becomes  very  strong,  and  Convection  B  becomes  very' 
weak  and  disappears  (Fig.5a)  due  to  the  high  growth 
rate  of  Convection  C  and  the  warm  SST.  The  certainty 
of  disappearance  of  Convection  B  makes  El  Nino  cycle 
relatively  regular.  The  strong  surface  easterlies  (from 
both  the  mean  winds  and  Convection  C)  switch  the 
OML  back  to  the  strong  surface  wind  forcing  regime; 
and  La  Nina  starts  again. 
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Fig.  5.  F.l  Nino  and  La  Nina  cycle. 


7.  SUMMARY 

This  study  shows  that  a  realistic  thermodynamics 
for  both  atmosphere  and  ocean  is  needed  in  the  ENSO 
theories.  As  the  ocean  surface  is  under  weak  wind  forc¬ 
ing,  the  time  rate  change  of  SST  perturbation  is  nega¬ 
tively  correlated  to  the  perturbation  of  the  ocean  mixed 
layer  (OML)  depth.  -  h'w  However,  as  the  ocean  surface 
Ls  under  strong  wind  forcing,  the  time  rate  change  of  SST 
perturbation  is  positively  correlated  to  the  perturbation 
of  the  OML  depth,  h’w.  Such  a  difference  leads  to  the 
generation  of  two  different  low  frequency  (interannual) 
modes  propagating  eastward  (strong  surface  wind  forc¬ 
ing)  or  westward  (weak  surface  wind  forcing).  A  new 
hypothetical  theory  about  La  Nina  El  Nino  cycle  is 
presented. 
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